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Q . We show that grains can be efficiently aligned by interacting with a subsonic 

j^. ! gaseous flow. The alignment arises from grains having irregularities that scatter 

atoms with different efficiency in the right and left directions. The grains tend to 
i — i align with long axes perpendicular to magnetic field, which corresponds to Davis- 

Greenstein predictions. Choosing conservative estimates, scattering efficiency of 
q impinging atoms and conservative "degree of helicity" , the alignment of helical 

£j ! grains is much more efficient than the Gold-type alignment processes. 
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00 ■ 1. Introduction 
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One of the most convenient ways to trace magnetic field is related to emission and ab- 
! sorption of polarized radiation by aligned dust. Therefore aligned dust is widely used for 
this purpose both in diffuse interstellar medium and molecular gas. Moreover, there is evi- 
.£h ! dence of aligned dust around young stellar objects and evolved stars as well as astrophysical 
environments. At the same time, the processes that aligned dust and their relation to mag- 
netic field are still the subject of debates and require further studies (see Lazarian 2007 for 
a review). While radiative torques (Dolginov & Mytrophanov 1976, Draine & Weingartner 
1996, 1997, Weingartner & Draine 2003, Cho & Lazarian 2005, Lazarian & Hoang 2007, 
Hoang & Lazarian 2007) are currently seen as the most promising candidate mechanism, the 
variety of astrophysical conditions may enables other mechanisms to dominate in particular 
environments. 

Mechanical alignment was pioneered by Gold (1951, 1952) and then quantified and 
elaborated by other researchers (e.g. Purcell 1969, Purcell & Spitzer 1971, Dolginov & 
Mytrophanov 1976, Lazarian 1994, 1997, Roberge et al. 1995). While the original mech- 
anism could deal with thermally rotating grains only, two modifications of the mechanism 
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introduced in Lazarian (1995) and elaborated later in Lazarian & Efroimsky (1996), Lazar- 
ian, Efroimsky & Ozik (1997), Efroimsky (1999) enabled the alignment of grains rotating 
at much higher rates. The latter were introduced to the field by Purcell (1979) (see also 
Lazarian & Draine 1999ab, where the limitations on the size for suprathermally rotating 
grains are discussed). 

The main shortcoming of the mechanical alignment processes was that they required 
supersonic gas-dust drift to get any appreciable degree of alignment (see Purcell 1969). 
Although later studies indicated that such drifts can be produced by ambipolar diffusion 
(Roberge & Hanany 1990, Roberge et al. 1995) or interactions of charged grains with MHD 
turbulence (Lazarian 1994, Lazarian & Yan 2002, Yan & Lazarian 2002, Yan, Lazarian & 
Draine 2004), the degree of alignment that is achievable for the Mach number drifts of the 
order of unity is insufficient to explain observations (see estimates in Lazarian 1997). 

The possibility of subsonic flows to mechanically align grains was mentioned in passing in 
Lazarian (2007) and Lazarian & Hoang (2007, henceforth LH07), but no relevant calculations 
were provided there. This possibility is related to irregular grains demonstrating helicity, 
i.e. ability of spin up in a regular way while interacting with a flow of particles. 



2. Toy model of a helical grain 

In LH07 radiative torques that arise from the interaction of photons with a grain were 
considered. Similar torques, however, should emerge when atoms bombard the surface of an 
irregular grain. 

In Figure 1 we provide a simple modef] of a helical grain. The grain consists of an 
oblate spheroid or an ellipsoid with a mirror attached to it at an angle 7r/4. The difference 
with LH07 is that we assume that the mirror is reflective to atoms rather than photons. 

For our purposes in the paper, the actual properties of a spheroid or an ellipsoid do not 
matter. Those geometric shapes are too symmetric and do not produce torques that can 
spin up a grain in a regular way, i.e. with angular velocity growing in proportion to time. 
The stochastic, i.e. with angular velocity growing as a square root of the time, spin up is 
the essence of the Gold-type alignment processes and it is ignored here, where we deal with 
subsonic flows. For the sake of simplicity, we assume that the Larmor precession of the grain 



1 Historically, simple models played an important role for developing models of mechanical alignment. 
The classical papers by Gold (1951, 1952) demonstrated the ability of grains to align using a thin stick as a 
model of a grain. 



- 3- 



in the external magnetic field is much faster than precession arising from the interaction of 
the ellipsoidal body with the gaseous flow. As a result, following arguments similar to those 
in LH07, one can disregard the effects of the gaseous flow on the ellipsoidal body altogether. 
Thus, the only torques to consider are those arising from the mirror. 

As in LH07, the the major role of the ellipsoidal body is to provide a steady rotation 
about its axis corresponding to the maximal moment of inertia. This ensures that the model 
grain is subject to regular torques while it interacts with the gaseous flow. 

The difference between the calculations in LH07 and those in the current paper stems 
from the fact that, while in LH07 the photons are coming as a beam from a single direction, 
for subsonic flow that we deal with, atoms hit the mirror from all directions characterized 
by the flux of atoms given by equation (j2J). This induces averaging of the torques that we 
implement below. 



3. Torques on a helical grain 

Consider a helical grain (see Fig. [I]) drifting across a gas chamber with a velocity 
along the e x axis. The velocity of an atom with respect to the grain is v — v d where v is the 
thermal velocity of the atom. The torque from the perfect reflecting of atoms on a surface 
area A of the mirror is given by 

dT = [r x AP]A/(s - s d )d 3 s, (1) 

where s = v/vth, vectors r = lif is the radius vector directed from the spheroid to the mirror, 
s d = Vd/vth and A is the surface area of the mirror, and 

f(s-s d ) = nv th (s-s d ).Ne- s \ for (s - s d ).N < 0, (2) 

is the flux of incoming atoms that can collide with the grain. Angular momentum element 
AP is antiparallel to the normal vector and given by 

AP = -2m H N\v - v d | cos 7, (3) 

where cos 7 = — (v — v^).N/|v — v d \. Here N are radius and normal vector which are 
functions of angles G,/3 describing the orientation of the grain in the lab system defined 
by ei,e2 and 63 (see Fig. [1]), as given in LH07@ Substituting equations (jSJ) and (j3J) into 
equation (JTJ we get 



2 N is given by equation (B6) in LH07, where n\ = —sin a and n 2 = cos a are components of N in the 
grain coordinate system. 
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dT = 2nm H v 2 h [(r x N]((s - s d ).N) 2 Ae- s2 d 3 s, (4) 

Taking into account the contribution from the reflection on the other surface, total torque 
becomes then 

dT = 2nm H v 2 h [rxN] / ((s-s d ).T<i) 2 Ae~ s2 d 3 s-2nm H v 2 h [rxN] / ((s-s d ).N) 2 ^ e - s2 (i 3 s, 

J(s-s d ).N<0 J(s-s d ).N>0 

(5) 

Integral (jSJ) can be analytically evaluated and the resulting torques are given by 
r= mH ™* hA [rxN]{\N 1 \N 1 (2s d e- s *+^erf 

(6) 

where Ni, N% and N$ are components of N in the lab coordinate system. 



Equation (J6j) can be written as 

m H nv 2 h Al 



-Or, (7) 



2 

where Qr is the vector torque efficiency consisting of components Q e i, Q e2 and Q e3 on axes 
ei,e 2 and e 3 , respectively. 

Two first components for subsonic and supersonic cases are shown in Figure [2j The 
third component produced by the mirror is found to be zero. It is clear that the component 
Qe3 is responsible for the precession (LH07). This component is non-zero for a spheroidal 
body of the grain. In fact, as we discussed above, we disregard this component, as it does 
not induce either alignment or spin up. 

In Figure [2] we show torque components for subsonic case M = 0.3 and a supersonic case 
M = 10. It can be seen that the component Q e2 has zero points at = 0, tt/2 and 7r in both 
subsonic and supersonic case. The component Q e \ is very symmetric in the latter. However, 
the torques are distorted due to the thermal collisions of gaseous atoms in the former case. 
The essential properties (i.e. symmetry of Q e \ and zero points of Q e2 ) of mechanical torques 
which are similar to those of radiative torques indicate that mechanical torques can align 
grains in the same ways as radiative torques do, i.e., grains tend to aligned with long axes 
perpendicular to magnetic fields. We study this problem in the section below. 



4. Grain alignment with respect to magnetic field B 



When the Larmor precession of a grain is larger than the precession induced by mechan- 
ical torques (see §11.6 in LH07), the alignment occurs with respect to magnetic field (either 
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|| or _L) which defines the axis of alignment. To study the grain alignment induced by me- 
chanical torques, similar to LH07, we solve the equations of motion for angular momentum 
J in time t: 

where T is given by equation ([6]) and t gas is the gaseous damping time. For simplicity, we 
assume a perfect internal alignment; i.e. the axis of major inertia a! is always parallel to 
J (DW97; LH07). We consider a model grain drifting across magnetic field by an angle 
ip = 70° acted upon by a subsonic flow with M = 0.3. Representing equation (jSJ) in the 
spherical coordinate system defined by J,£,(j) where J is the magnitude of J, £ is the angle 
between J and B, and is the Larmor precession angle , we get three equations for J, £ and 
4> (see Sec. 7 in LH07). Averaging over <fi from to 2tt, we obtain only equations for J and 
£. We construct trajectory map based on the obtained solution J(t),£(t). 

Figure [3] shows the obtained trajectory map with initial condition J/Jth = 200 and the 
angle £ generated from a uniform distribution. We see that grain phase trajectories that 
start at different angles £ converge on either high- J attractor points A and B corresponding 
to perfect alignment, or a low-J attractor point C. The latter point formally corresponds 
to J = 0. However, as it was explained in LH07 and shown by numerical simulations in 
Hoang & Lazarian (2007), this is an artifact of our ignoring thermal fluctuations within the 
grain. If those are taken into account the low-J attractor points correspond to the angular 
momentum of the order of the thermal angular momentum corresponding to the temperature 
of the grain body (see also Weingartner & Draine 2003, Hoang & Lazarian 2007). 

What should be noted is that the attractor points for M = 0.3 correspond to the 
alignment of &i parallel to B; i.e. the grain gets aligned with long axes perpendicular to 
magnetic field. This is similar to what the Davis- Greenstein (1951) mechanism of param- 
agnetic relaxation predicts, but the process we discuss here does not invoke paramagnetic 
relaxation. 

Gaseous bombardment randomizes grain phase trajectories. However, rather than trac- 
ing phase trajectories, it is practically convenient to have a criterion for the efficiency of 
grain alignment. The maximal rotational velocity of grains induced by torques can be used 
for this purpose. A study in Hoang & Lazarian (2007) which solves the Lagevin equation in 
order to account for gaseous bombardment shows that this velocity which us a function of 
ip should be approximately 3 times larger than the thermal rotational velocity of the grain. 

For practical calculations of the rotational velocity we assume that 10% of the atoms 
impinging on the grain mirror are reflected. Figure H] shows the angular velocity of a grain 
as a function of the Mach number of the flow. The horizontal line corresponds to J ma x — 
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3IlU>therm,,gas for ij} — 0°. 

It is clear from looking at Figure H] that for our model, helical grains get aligned for 
relatively low velocities of gas-grain drift. All the earlier mechanisms of mechanical alignment 
are inefficient for such low subsonic velocities. Indeed, drifting velocities of the order of 0.1 
of sound speed should be quite common in the ISM (see Yan, Lazarian & Draine 2004), so 
potentially the mechanism is widely applicable. However, it requires further studies. 

One may wonder whether the condition of perfect reflection is absolutely necessary. 
It can be shown that the torques are changed by a factor of unity if impinging atoms are 
absorbed on the grain surface, thermalized there, and emitted from the point that they hit 
the surface. Moreover, for regular torques to act on a helical grain, it is sufficient to have a 
correlation of the place that the atom impinges on the grain surface and is evaporated from 
it. It is only when there is no correlation at all, that our model grain does not experience 
regular torques. Therefore we believe that we do not overestimate the torques with our 
assumption of 10% reflection. In general, to characterize the interaction of grain surface 
with the gaseous flow, one can introduce a "reflection efficiency factor" E. 

A different issue is the degree of helicity of an irregular grain. Naturally, most of the 
grains do not have facets, but have numerous irregularities. An idealized model grain (see 
Fig. 1), for which the damping arising from the gas interactions with the ellipsoidal body 
tends to zero, has the maximal possible value of the torques Q opt (a, M) for a given grain 
size a and the Mach number M and its "helicity reduction factor" D — 1. This factor, in 
a general case, can be characterized by D = J? , M )^ . Further numerical studies should 
clarify the value of D for actual irregular grains. In Figure H] we show results for different 
ED factors. 



5. Discussion and Summary 

In our Letter we use the same toy model of a helical grain as in LH07. In LH07 we 
proved the validity of the model (which uses a geometric optics approximation) comparing the 
functional dependencies of model's torques and those numerically calculated for the actual 
irregular grains. In comparison, the justification of the model for the gaseous bombardment 
is self-evident. However, the combination ED presents a combined uncertainty factor, which 
characterizes both the uncertainties in grain-gas interactions and grain shape. Testing of 
the former requires laboratory studies, while the latter can be established via the numerical 
research. 

What is the relative role of traditional mechanisms of mechanical alignment? Those 
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mechanisms require supersonic motions to be efficient. It is only mechanical alignment of 
helical grains that allows alignment for the subsonic motions. Thus we may claim that 
mechanical alignment of helical grains constitutes a new class of mechanical alignment pro- 
cesses. In fact, the difference between the mechanical alignment of helical grains and the 
Gold mechanism of grain alignment is similar to the difference between the radiative torque 
alignment and the Harwit (1970) process. The latter appeals to stochastic torques arising 
from absorption of photons from a light beam and was shown by Purcell & Spitzer (1971) 
to be inefficient for most interstellar situations. 

Similarly, as the radiative torque alignment is more efficient than the Harwit (1970) 
alignment, we believe that the mechanical of helical grains is more efficient than the Gold 
processes. Therefore, at supersonic velocities we expect the mechanical alignment of the 
actual irregular grains to be governed by their helicity, rather than the degree of their 
oblateness of prolateness, which is the case for the Gold mechanisms. The consequence of this 
is that the mechanical alignment of grains will happen with long grain axes perpendicular 
to magnetic field irrespective on the gas-grain velocities. This is in contrast to the Gold 
alignment, for which the change of alignment happens at the Van Vleck angle a = 54°. 7. 

In Hoang & Lazarian (2007), the effects of gaseous bombardment and uncompensated 
Purcell's torques (Purcell 1979), e.g. H2 torques, was studied in connection with the align- 
ment of helical grains subjected to radiative torques. That study shows that, in the case 
when both high-J and low-J attractor points coexist, gaseous bombardment transfers grains 
from low-J to high-J thus increasing^], counter-intuitively, the degree of alignment. This is 
surely true only if J m ax(ip) is sufficiently high, e.g. higher than 3Iiu t h e rm,gas- The Purcell 
torques, do not change J in the low-J attractor points because grains flip fast at those points 
causing thermal trapping, as discussed in Lazarian & Draine (1999a). They, however, can 
increase the values of J at high-J attractor points. For most situations, this does not affect 
the degree of alignment. 

Our considerations above deal with ordinary paramagnetic grains. For those, the in- 
fluence of paramagnetic torques is mostly negligible for the typical ISM conditions. The 
situation changes if grains are superparamagnetic (see Jones & Spitzer 1967). In this sit- 
uation one can expect always to have alignment with high-J attractor points, thus having 
degree of alignment close to 100%. The assumption of the presence of superparamagnetic 
inclusions is an extra assumption, however. 

As we have stressed above, the mechanical alignment we discuss here and the radia- 



3 The increase mostly stems from the higher degree of internal alignment in high-J attractor points 
compared to the low-J ones. 



- 8- 



tive torque alignment are different incarnations of the alignment of helical grains. While 
the radiative torques have attracted a lot of attention recently, the mechanical alignment 
of helical grains has only be mentioned in a couple of publications (Lazarian 1995, 2007, 
Lazarian et al. 1997, LH07). The relative role of the mechanisms depends on the yet un- 
certain factor ED for the mechanical processes. The observations tends to be in agreement 
with the radiative torque predictions (see Lazarian 2007 and ref. therein). It is encouraging 
that both mechanisms predict the alignment with long grain axes perpendicular to magnetic 
field. Therefore one may hope that the subsonic alignment of helical grains can reveal via 
polarimetry the magnetic fields in the situations when the radiative torques fail. The variety 
of astrophysical circumstances ensures that there are situations when grains are aligned by 
mechanical subsonic flows. Note that aligned grains were reported not only for the ISM, 
but also comets^ and circumstellar regions. They are also likely to be present in the disks 
around evolved and young stars. 

The goal of this Letter above is to attract the attention of the community to the pos- 
sibility of mechanical alignment of helical dust grains. While further studies of the process 
are necessary, the following points can be made at this moment: 

1. Irregular grains are, in general, are expected to exhibit helicity when they interact 
with gaseous flows. 

2. The mechanical alignment of helical grains is efficient even when the flow of gas is 
subsonic. This alignment is likely to dominate Gold-type processes. 

3. The mechanical alignment of helical grains aligns grains with long axes perpendicular 
to magnetic field, thus increasing the chances of magnetic field tracing when other alignment 
mechanisms fail. 

Acknowledgments. The work was supported by NSF grant AST 0507164 and the 
NSF Center for Magnetic Self-Organization in Laboratory and Astrophysical Plasmas. 
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Fig. 1. — A toy model of helical grain consisting of a reflecting spheroid and a mirror attached 
to it end in gas flows. The distance between the mirror and the spheroid, l\ is assumed to 
be much larger the the mirror size li- 



1.5 



Cr 




-1 .0 



-1 .5 I- . i .... i .... i .... - 

-1.0 -0.5 0.0 0.5 1.0 

COS0 

Fig. 2. — The normalized torques Q e i(©) and Q e 2(0) averaged over the rotation are shown 
for different Mach numbers M of the relative grain-gas flow. For large M = 10 the torque 
Q e i(0) is symmetric and Q e 2(©) has three zeros at O = 0,7r/2 and 7T. For subsonic flows 
M = 0.3, the torques are modified due to the effect of thermal velocity, but zeros points of 
Q e2 (6) are unchanged. 
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Fig. 3. — Phase trajectory of the grain of li = IO/2 = 0.1/xm and ED = 1CT 3 , in diffuse gas 
with n = 30cm -3 , temperature T = lOOfT corresponding to the March number M = 0.3 for 
an angle ip between and magnetic field B equal 70° has attractor points A and B of high 
angular momentum, C with J = . In the presence of thermal grain wobbling the point C 
becomes a low- J attractor point with J ~ Jtherm,dust- 
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Fig. 4. — Ratio of the angular momentum induced by grain drifting to the thermal value of 
the model grain as a function of March number for different "helicity reduction factors" D. 
The horizontal line corresponds to J max (V ; ) = 3Ii^therm,gas which presents the threshold for 
grain alignment obtained in Hoang & Lazarian (2007). 'The results are obtained for a gaseous 
flow moving parallel to the magnetic field. For grains with velocity mostly perpendicular 
to the magnetic field (see Yan & Lazarian 2003), the angular momentum J can be further 
decreased by a factor 10-100 



